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Current errors on \V u b\ are dominated by model dependence. For in- 
clusive decays, the model dependence comes from the Fermi motion of the 
b quark. By combining the endpoint photon and lepton spectra from the 
inclusive decays B — > X s 7 and B — > X u t 77, it is possible to remove this 
model dependence. We show how to combine these rates including the 
resummation of the endpoint logs at next to leading order. The theoreti- 
cal errors on \V U }\ on the order of 10% are possible. We also give a brief 
discussion on comparing different extractions. 
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1 Introduction 



The Cabibbo-Kobayashi-Maskawa matrix element V u b is very important for un- 
derstanding CP violation in the Standard Model. An accurate measurement of \V u b\ 
puts strong constraints on the Unitarity Triangle. Unfortunately, V u j, is vary hard 
to measure. Current measurements have errors that are dominated by model depen- 
dence. Some of the best extractions so far have come from exclusive decays, such as 
B — ► tt£v or B — > ptP. The problem with exclusive decays is the strong hadronic dy- 
namics can not be calculated, and we have to resort to models, light-cone sum rules, 
or lattice QCD calculations to obtain the form factors [Q. At the present time, all 
these methods give around 20% errors. A recent measurement from CLEO || using 
B -> plv gives \V ub \ = [3.25 ± 0.14(stat.)^(syst.) ± 0.55(model)] x KT 3 . In the 
future, the lattice will give accurate predictions for the form factors, but until then, 
a measurement of 20% is probably the best we can hope for from exclusive decays. 

In some ways, inclusive decays should provide a straightforward means to measure 
\V u b\. All we need to do is measure the total rate b — *> u£v, which is proportional to 
\V u t,\ 2 and is known to order a 2 s @. If we could measure the total rate, we would not 
have to worry about quark-hadron duality violations, thus a very accurate measure- 
ment would be possible. 

Unfortunately, there is a very large background from h — > c decays, which is about 
100 times more abundant than b —>■ u decays. To remove this large background, 
kinematic cuts must be made. Three basic cuts are discussed in the literature, each 
having its own advantages: a cut on the electron energy spectrum, a cut on the 
hadronic invariant mass spectrum, and a cut on the leptonic invariant mass spectrum. 
For now we will concentrate on the electron energy spectrum, and return to the other 
cuts later. 

Since the u quark is much lighter than the c quark, the electron energy spectrum 
for b — > u decays extends past the endpoint for b — > c decays, see Fig. [I]. Thus 
it is possible to remove the charm quark background by cutting above the b — > c 
endpoint. All that is necessary is a theoretical prediction for the integrated rate 
above the cut. Unfortunately, putting a cut near the endpoint introduces a new 
small mass scale, AE ~ 300 MeV, which introduces large perturbative [log(rrib/AE)] 
and non-perturbative (A/AE) corrections. Therefore both the perturbative and non- 
perturbative series must be resummed for the rate to be trustworthy. 

The calculation of the rate begins with the effective Hamiltonian Q 



obtained by integrating out the t quark and W bosons. The differential decay distri- 
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Figure 1: Electron spectrum for semi-leptonic b decay rates to c and u quarks. The rate 
for b — > ueV has been multiplied by a factor of 10. The region in AE can only have b — ► w 
decays, and thus is useful for extracting |V^6|. 



bution can then be written as the product of leptonic and hadronic tensors 

dT oc \V ub \ 2 L a ?W a(3 . (2) 

Using the Optical Theorem, the hadronic tensor W a p can be related to the imaginary 
part of the time ordered product of currents 

W aP = —ImT af3 , (3) 

7T 

T aP = -^-Jd"xe-^{B\T(4(x)M0))\B). (4) 

The time ordered product can be calculated by expanding in an Operator Product 
Expansion (OPE). The Wilson coefficients can be calculated (over most of phase 
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space) in perturbation theory while higher dimensional operators in the OPE 

are suppressed (over most of phase space) by powers of l/m b . Thus, the leading term 
in the OPE gives free b quark decay. The first corrections enter at order 1/mf, and 
are proportional to the Heavy Quark Effective Theory parameters 
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The problems begin as the energy of the lepton approaches the endpoint. Defining 
x = 2Ei/mi ) to be the rescaled lepton energy, the higher dimensional operators in the 
OPE are actually suppressed by 

A 
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Higher dimensional operators in the expansion are no longer suppressed. In other 
words, the expansion is becoming singular as we approach the endpoint. 

The breakdown in the OPE can be seen in the expression for the rate at order 
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by the appearance of singular functions at the endpoint. 

To handle the breakdown of the non-perturbative series, the leading singular terms 
must be resummed. These corrections resum into a non-perturbative structure func- 
tion, f(k + ) |J. The differential rate is now a convolution of f(k + ) with the partonic 
rate Hi 

' /J f dk + f(k + )^(mt), (9) 



where ml = m^ + k^ 
function occurs for b 
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The structure function is universal, meaning that the same 
u£i7 and b — > 37 decays. Being a non-perturbative function, 
f(k + ) is not known; we do know the first few moments of f(k + ), however. Thus, 
to handle the endpoint region, some model for f(k + ) must be introduced. We could 
in principle extract the structure function from b — > sj decays and then apply it to 
b — > uiu, but this is difficult because of the way f(k + ) enters the rate (||). Instead, 
we will skip the step of extracting the structure function and directly use the b — > S7 
rate in the b — > uiv rate. But first we need to discuss the perturbative corrections. 



Near the endpoint, the perturbative correction to the rate looks like 1 1 
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As x — ► 1, the logs become large and the perturbative series breaks down. To trust 
the prediction, the logs need to be resummed. There are similarly large logs in the 
rate for b —>■ sj, so the logs must be resummed there, too [ |TT|1 . 

It is possible to resum the series using Infrared Factorization, which is also used 
for DIS, Drell-Yan, etc. The idea is that in the endpoint region, the light quark is 
shot out with large energy, but with small invariant mass. This quark produces a 
jet of particles through collinear radiation. While the constituents of the jet can talk 
to each other (and the original b quark) through soft gluons, hard gluon exchange is 
disallowed. The soft radiation cannot tell if the jet was initiated by a u quark or an 
s quark, thus it will be the same for b — > u and b — > s decays. 



Mathematically, there is a separation of momentum regions into [12 



Hard (H) 
Jet (J) 
Soft (5) 



k+~k_~k t = 0(m b ), (11) 
k + = 0[m b (l - x)], fc_ = 0(m b ), kt = 0(m b Vl - x), (12) 
k + ~k_~k t = 0[m b (l - x)}. (13) 



By introducing a factorization scale \x to keep these regions separated, we can write 
the rate in factorized form as 



dT 

dx 



dzS(z,fj)J(z,n)H(n). (14) 



The soft function S(z,fM) is the same for b — > u and b — » s, while J(z,fM) and H(/i) 
depend on the process. 

The rate completely factorizes after taking moments, 



Mjj = I dxx n - 1 ^ f — = S N J] r H] T , (15) 
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J dx v S N J s *{x v )H%{x v ), (16) 



where and are the moments of the b — > and b — > u£u rates, respectively. 
The soft function contains perturbative and non-perturbative pieces 
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Inctn, (17) 



where are the moments of the structure function introduced earlier. Thus we can 
write the moments fllSD and (|16]) as 

Ml = f N a N Jj,Hl, (18) 

M§ = J dx v f N a N J$(x v )H$(x v ). (19) 
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All the large logarithms are contained in the combination (JnJn- The only fact we 
need about the perturbative resummation is that after resumming, including next- 
to-leading logarithms, there is the relation |y| 

a N J$ = a N J]fexjp[g sl (a s logiV)], (20) 

where exp[g s i(a s log N)] is a known function. 

We can now combine the above results. Substituting first (|20| ) into (|1^), and then 



(P~8| ) into the result, we get 

M$ = / dx v f N a N J^{x v )H s ^{x v ) 



dx v f N a N Jjj exp [g si (a s log N)]H S N (x 
M 



f 1VL 

= J dx v -j^e W [g sl {a s \ogN)\H s ^x v ). (21) 

Note that the dependence on the unknown structure function has been eliminated. 

We can go back to x-space by taking an inverse Mellin transform. The left-hand 
side of ( [Zip is just the semi-leptonic rate. The right-hand side is a convolution of the 



b — > S7 rate with a known function. Rearranging, we can write this as JTJ] 

\V ub \ 2 JTjb^uiu) 
\V t * s V tb \ 2 JfdTydx~<*K(x~<,a s 



[22) 



So in words, what we have done is written \ V u b\ 2 as the ratio of the b —>■ ueV rate over 
a convolution of the b — > 37 rate with a known function. 

What are the uncertainties? First, there are higher order corrections that we 
neglected, which enter at the order of A/m^, a s (l — x) and (1 — x) 3 . For the value 
of the electron energy cut, x cut ~ 0.87, these corrections should all be less that 10%. 
Of course, we are estimating the size of the higher order corrections, since they have 
not been calculated. They may be larger or smaller by a factor of 2 or 3. Without 
calculating the corrections directly, it is not possible to know. We will come back to 
this qualification shortly. 

Second, there are the violations of quark-hadron duality. These violations are 
hard to quantify, but they should be small if we are not dominated by a just a few 
resonances; the more final states, the smaller the duality violations. In the region 
that we are interested in, it does not appear that we are dominated by resonances, so 
neglecting them should be okay. It would be better if we could have a larger number 
of the decay products. 

This is possible if we cut on different kinematic variables. The other variables 
discussed in the literature are the hadronic invariant mass ||15|| , and the lepton in- 



variant mass |L6] . The hadronic invariant mass spectrum also has dependence on the 
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structure function |l5|Jl7| , which introduces model dependence. However, by using a 
method analogous to the one described above for the electron spectrum, the depen- 
dence on the structure function can be eliminated [[H|,[n|. The errors from higher 
order corrections are similar to the electron spectrum and should be around 10%. 
The main advantage of the hadronic invariant mass is that after a cut to remove the 
charm background, between 40% and 80% of the possible final states will be included. 
This is much larger than for the electron spectrum, which includes about 10% of the 
possible final states. Thus the quark-hadron duality violations should be negligible. 

The leptonic invariant mass cut has different advantages [16,2(J. Here the struc- 
ture function is not important, so we do not need to do anything to remove this model 
dependence. Higher order non-perturbative corrections are on the order of (A/m c ) 3 , 
which leads to an error again of around 10%. This disadvantage for this cut is the 
fraction of final states included after the cut is around 20%, so the quark-hadron du- 
ality errors may be an issue. There is also some question about how good a resolution 
can be obtained on the lepton invariant mass, which is the only immediate problem 
for this method. 

All three of the above methods should have a theoretical uncertainty (modulo 
quark-hadron duality violations) of around 10%. Again, these are estimates of higher 
order corrections. The actual errors may be bigger or smaller. Also, the duality 
violations could enter in different ways for each measurement. To really trust any 
extraction of \V u b\, we should measure it as many ways as possible, and only after (or 
if) there is a convergence of the results should we trust the extracted value. 
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